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We  consider  the  lattice  gas  model  of  the  0/W(110)  chemisorption  system. 
We  calculate  the  transition  temperature  T^  at  one  half  coverage  by 
means  of  a  new  renormalization  group  treatment.  Good  agreement  10%) 

'  t 

with  Tj,  obtained  via  Monte  Carlo  studies  1$  found  over  a  wide  range 
of  parameter  values.  We  explore  the  remarkable  independence  of  Tj.  from 

the  nearest  neighbor  adatom-adatom  interaction.  We  also  give  several 
arguments  indicating  that  the  transition  at  this  coverage  is  in  the  simple 
Ising  universality  class. 
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1.  Introduction 


The  Interaction  between  adoorbed  atoms  Is  a  fiindamental  property 
of  adsorption  systems.  The  understanding  of  many  Important  surface  phenomena 
such  as  corrosion,  oxidation,  hetergeneous  catalysis  and  crystal  growth  will 
be  enhanced  by  a  knowledge  of  these  forces.  In  order  to  fully  understand 
adatom-adatom  interactions  It  Is  necessary  to  start  with  a  study  of  adatom 
behavior  in  submonolayer  systems.  This  allows  the  isolation  of  the  effect 
from  the  many  other  phenomena  occuring  in  adsorption  systems. 

Much  progress  has  been  made  in  the  ab  initio  calculation  of 

adatom-adatom  (AA)  interactions,  especially  for  physi^sorbed  layers. 

Here  we  consider  the  strong  chemisorption  system  0/W(110).  In 

this  kind  of  system  indirect  AA  interactions,  mediated  by  the  substrate 

conduction  electrons,  are  important.  A  general  theory  of  these  forces  existj/ 

but  in  order  to  determine  AA  couplings  for  any  specific  case  it  is  necessary  to 

proceed  by  a  more  indirect  route.  In  the  0/W(110)  system  models  for  the  AA 

interactions  have  been  derived  by  making  use  of  statistical  mechanics.  In  the 

submonola^er  regime,  for  T  ^  300°K,  the  0  overlayer  is  in  thermodynamic 
H- 

equilibrium.  Hence 'its  thermodynamic  phase  is  determined  by  the 
A-A  interactions.  The  phase  diagram  of  this  system  can  be  measured  with  a 
combination  of  LEED  (Low  Energy  Electron  Diffraction)  and  Auger  techniques.  Using 
statistical  mechanics,  one  can  then  compare  the  experimer  base  diagram 
with  one  calculated  from  an  assumed  set  of  AA  coupling  cc  sii This  can  be 
done  for  various  parameter  values,  and  the  best  AA  interaction  energies  may 


Acc«ssl/n 


be  determined.  This  program  has  been  carried  out  for  0/W(110)  using  the  Monte 


»TIS  >  Carlo  method.  This  in  a  very  powerful  numerical  technique,  but  its  repeated  use 
DTIC  T.5B 

Unanaour^'  requires  nun-trivial  amounts  of  computer  time.  In  this  work  we  show  that 
Justifies 

■  Comparable  results  may  be  obtained  more  efficiently  with  a  renormalization  group 

By - - 

Distrlb’.t  method.  Thin  type  of  calculation  was  originally  applied  by  Wilson  to  the  problem 
Availal ■ 

of  determining  thermodynamic  behavior  near  critical  points.  It  has  more  recently 

Dist  ?:p(.  cl.i.  I 
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been  recognized  as  also  providing  a  powerful  tool  for  approximating  the  global 
free  energy  and  phase  behavior. 

In  this  work  we  apply  an  existing  position  Bjiace  renormalization  group 

(PSRG)  technique,  the  modified  Kadanoff  variational  method  (MKVM)  to  the 

problem  of  0/W(110).  The  model  used  for  this  adsorption  system  Is  a  lattice 

gas  with  Interactions  through  4th  nearest  neighbor.  Thus  our  work  represents 

7 

a  significant  extension  of  existing  PSRG  methods.  For  the  case  of  one  half 
coverage  we  find  excellent  agreement  with  previously  reported  Monte  Carlo 
results  over  a  wide  range  of  parameter  values. 

In  section  2  we  introduce  the  lattice  gas  model  for  our  adsorption 
system  and  describe  a  certain  useful  approximation  we  have  developed  for  It. 

Section  3  describes  the  MKVM  as  applied  to  this  system.  In  section  4  we 
present  our  results  and  compare  them  with  those  obtained  by  the  Monte  Carlo 
method.  Section  5  we  discuss  the  critical  properties  of  the  model.  We  point  out 
that  there  is  considerable  evidence  that  the  order-disorder  transition  at  one-half 
coverage  is  Ising  like.  In  section  6  we  mention  some  problems  for  future  research. 


2 .  Definition  of  Model  and  the  Approximation 

The  adsorption  sites  for  0  on  the  W(llO)  have  the  geometry 
illustrated  in  Fig.  1.  Note  that  this  two-dimensional  lattice  may  be 
regarded  as  two  interpenetrating  rectangular  sublattices  R  and  R'. 

The  Hamiltonian  for  this  system  may  be  written  as 

It  -  '5-  r  V  %  T.  i'  ^  2,  r-  +  c-  j  ?.  r, 


Tlie  occupation  number  ^  =  1  (0)  for  an  occupied  (empty)  adsorption  site. 

The  /\A  energies  are  defined  with  respect  to  the  (isolated)  single  adatom 

adsorption  energy  Tliis  is  of  the  order  of  several  eV  and  plays  no  role 

at  the  tL  rnperatures  of  interest  here  since  kT t-  .  ^  nn'>  denotes  a  sum 

o 


over  all  distinct  nearest  neighbor  adsorption  site  pairs,  H  over  (horizontal) 
second  neighbor  pairs,  v  over  (vertical)  third  neighbor  pairs,  and  D  over  fourth 


neighbor  pairs  In  the  diagonal  directions.  The  grand  partition  function  Is 


'I- 


«.  ,  -p‘U(C.-TJ 

U,W;T)r  51  i  2.'  e  ,• 


(2) 


^  represents  a  sum  over  all  possible  occupation  numbers  satisfying 


(3) 


and  removing  the  prime  removes  this  restriction.  N  is  the  number  of  adsorption 


sites.  Eq  (2)  may  be  rewritten  using  spin  variables 

l-\  ^(r,  • "  57/ ) 


(4) 


where 


K=  -I  'h?:/2  4-  (^/V  6^4  Z6^)l 

-  L  l^i/2  4-  42  6i-»16j  +‘/6^  )]2<r 


(5) 


Here  ~  2*?^  ^  so  =4-1  (spin  up)  corresponds  to'T^  =  +  I  (filled) 


site  and  (j7  =  -  1  (spin  down)  to  =  0  (empty  site).  Now  Eq(5)  is  an  even 


function  of  (Q  when 


(6) 


Hence  imposing  Eq(6)  implies  that  the  coverage  C'  =  Na/N  =  1/2. 

Now  as  pointed  out  above,  the  adsorption  sites  for  0/W(110)  may  be  regarded 
as  two  interpenetrating  sublattices  R  and  R' ,  an  shown  in  Fig.  1.  Hence 
Eq(5)  may  be  written  as 
11 


(7) 


where 


\/Af(^)  -- '  f'/q"K5b. i  /t^) A'', fr.ifyTj 


(8) 


and  IP  (j'  ij",  are  spin  variables  for  one  cell  in  the  rectangular  lattice 
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R  as  shown  in  Fig.  2.  The  sums  on  R  and  R'  in  Eq(7)  are  over  distinct 
4-Bpln  sublattice  cells. 

Note  that  Is  the  only  A~A  Interaction  ta'na  coupling  sublattices 

R  and  R'  in  Eq(7).  Even  though  this  energy  is  as  large  or  larger  than 
the  further  neighbor  coupling  terras,  it  turns  out  to  have  very  little  effect 
on  the  value  of  the  (order-disorder)  transition  temperature  at  coverage 
6  =  J5.  This  was  already  noted,  for  a  certain  range  of  values,  by  Williams 
et  al  in  their  Monte  Carlo  study  of  this  model.  We  examine  this  indepence 
further  and  use  it  In  our  PSRG  treatment  below. 


Now  for  0/W(110)  the  ordered  (low-temperature)  state  at  ■  *5  is  (2x1). 

This  ordering  also  occurs  at  low  coverages  (although  the  transition  temperature 

r 

drops  abruptly  for  0.3).  These  facts  imply  that  one  may  feke 

<0  (9) 

(In  fact  one  could  also  let  62  and  slightly  negative). 

Consider  a  perfectly  ordered  (2x1)  state  at  ^  =  *5.  Moving  an  adatom  to 

the  most  favorable  disordered  site  raises  its  energy  by  0'>'' 

^  ,■  hr 

^quantity  independent  of  t-^.  This  implies  that  T^  depends  on  weakly  at 

best,  as  is  verified  in  Ref.  Z.  for  large  values  of  6-^. 

When  is  small  we  may  also  use  an  argument  following  Kadanoff  and 

Wegner.  Eq(9)  implies  that  the  sublattices  R  and  R'  Will  be  antif erromagnetlcally 

ordered  at  low  temperatures.  Hence  the  Hamiltonian  is  of  the  form  of  two 

interpenetrating  antiferromagnetic  square  lattices  coupled  f erromagnetlcally 

-  L  2  <r  T-  2 'r//r  '  -  i<,  2.  <r  r '  (lo) 

where  k,  kj^'^  0.  In  Eq(lO),  Q~^  are  spins  on  sublattice  R,  on  sublattice  R' , 
and  the  three  sums  are  over  nearest  neiglibor  pairs  both  on  R,  both  on  R'  and 

one  on  each,  respectively.  Now  for  =  0  we  have  two  independent  lattices 

with  an  antiferromagnetic  transition  at  coupling  K^.  The  correlation  function 
associated  with  the  operator  connecting  the  two  lattices  is 

extn)  ^  <u'.  r. ‘  cr  ,  r/  > 

J  *•  J  I J ,  r  / 


(11) 
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where  i  and  j  are  nearest  neighbors  on  R  and  R'.  For  =  0,  at  the 


critical  point  one  has 


-  ~  -k. 


To  evaluate  ^ ,  we  must  examine  the  large  r  behavior  of  the  ferromagnetic 
correlation  function  at  the  antlferromagnetlc  critical  point.  By  up-down 
spin  symmetry,  this  Is  the  same  as  the  antlferromagnetlc  (staggered) 
correlation  function  at  the  ferromagnetic  critical  point.  This  falls  off 
more  rapidly  with  r  than  the  ferromagnetic  correlation  function,  for  which 
{(ivj  =.  *?/*/  (in  two  dimensions).  Hence  we  must  have  x  >9/4.  Now  let 

Then  the  antlferromagnetlc  correlation  function  may  be  written  approximately 
as 

where  a  is  a  lattice  spacing.  For  v this  reduces  to 

Thus  X  =  17/^.  Now  y=d-x=2-  17/4  =  -  2  1/4  ^  0.  Hence  the 

operator  coupling  the  two  lattices  (the  third  term  in  Eq(7))  is  irrelevant 

at  the  fixed  point.  Tliis  ineans  that  for  small  it  should  have  a  vanishing 

affect  on  the  renormalization  group  flown  and  hence  not  affect  T^  . 

Since  thin  argument  depends  only  on  the  long  distance  behavior  of  certain 

correlation  functions  and  the  types  of  coupling  present.  It  Is  valid 

for  many  Hamiltonians.  Tn  particular,  Eq(7)  for  £r.=0  has  an  Islng  type  critical 

6.!u«  l.-'A. 

point.  Hence  ">]  =  1/4  and  the  argument  remains  valid.  We  use  it  below  to  justify 
treating  the  6^^  term  In  nq(7)  approximately. 

ij 

Actually  the  weak  cKpindcnce  of  T^  on  extends  over  a  wide  range  of 
(r  ^  values.  This  is  llluntr.ated  In  columns  4  and  5  of  Table  I  where  we 
list  T^  for  the  paraneter  values  in  Columns  1,  2  and  3  an  determined  by  Monte 
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Z.  w 

Carlo  calculations  (Col.  4)  and  for  =  0  as  given  by  the  free  fermion 

ej  ^ 

approximation  of  Fan  and  Wu  (col. 5).  The  uncertainty  in  col.  4  isilO  K  and 

10 

the  col.  5  values  are  about  32  too  large  for  the  {Parameters  shown. 


3.  Method  of  Calculation 

6 

In  this  section  we  describe  the  MKVM  technique  and  how  it  is 
applied  to  our  model  for  0/W(110).  This  is  defined  by  Eq(7) .  First  we 
consider  the  case  of  =  0,  so  that  the  two  sublattices  (see  Fig  1  and 
Section  2)  are  decoupled.  Later  we  include  via  a  preliminary  transformation. 
For  =  0  the  remaining  terms  in  Eq(7)  are  completely  Independent 


Eq(4)  reduces  to 


t  e 


Thus  the  free  energy  is  given  by  .  « 

- 


(17) 


We  calculate  ^nd  its  derivatives  using  the  modified  Kadanoff  variational 

method  (MK'.T'l)  •  With  the  signs  of  the  AA  energies  6-,^  defined  in  Eq(9),  an 
antiferromagnetic  coupling  is  favored  on  each  of  the  sublattices  R  and  R'. 

To  obtain  a  viable  cell  spin  definition  in  the  MKVM  (especially  at  low 
temperatures),  wo  must  redefine  spin  variables  so  that  ferromagnetic  ordering 
is  favored.  For  each  cell  in  sublattice  R  (and  likewise  for  R')  we  let 


are  Illustrated  in  Fig.  2.  Tlie  two  different  kinds  of  spin  redefinition  are 


performed  alternately,  so  that  the  configuration  sum  in  Eq(17)  is  equivalent 
to  a  conf  iguration  s\im  in  the  original  variables  with  ,  t  x  A' replaced 
by  ^  w  respectively.  Then  we  apply  the  MKVM  to  Eq(17) 

«lth  f-,,  replaced  by  -  i,  (-4  •  Thus  becomes 
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(18) 


(Note  also  that  this  transformation  changes  the  term  coupling  to  ()^  (0^  f  (Tj  4  41) 
into  a  "staggered  field"  term.  The  coefficient  of  this  term  vanishes  for 
Q h,  however.) 

6  , 

Now  we  review  the  MKVM  technique.  The  transformed  cell  potential  V 

for  the  new  lattice  with  doubled  lattice  spacing  is  related  to  the  original 
■<le 

cell  potential  V  by  RG  transformation  equation: 


where 


^vVyu)  ^  f 


(19) 


U.  (.(P,  <>•  (T^  )  -  1 2  Co5,t\  P  ( +  •j'l  i  'J)/  ) 

Z  =  no.  of  nearest  neighbors  =  4,  and  ^  spin  variables 

in  the  new  lattice  with  variational  parameter  p  is 

determined  by  minimizing  the  single  cell  free  energy.  This  leads  to  the 
nonlinear  equation 


2. 

-/Wy 


■  P'* 


.)V 


V(>A) 


V(ju) 


-  o 


(20) 


which  is  solved  at  each  iterative  step  to  determine  p. 

To  begin  with,  we  use  V  of  Eq. (18)  in  the  right  hand  side  Eq(19)  with  p 
given  by  Eq(20),  and  calculate  the  transformed  cell  potential  V  from  Eq(19) . 

Thin  constitutes  the  first  step  of  the  RG  transformation.  Note  that  this  step  in 
general  given  rise  to  a  non-zero  four  spin  coupling  term  (  aLjU,,^Uj^U^  Uy  ) 
well  as  the  terms  of  Eq(8).  In  addition  there  will  be  an  additive  constant 
term  in  V'  ^0.  We  then  use  V'^*)  thus  obtained  as  input  in  the  name  procedure 
to  calculate  the  transformed  potential  V"(i^.  This  RG  transformation  is 
iterated  further  so  that  a  series  of  cell  potentials  V'^),  V"^ 
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and  corresponding  variation  parameters  . P^*^  »  ••  obtained 

In  this  step  by  step  RG  transformation,  the  variational  parameters  p  '  and 
coupling  constants  contained  In  V'  '  '  (le  the  coefficients  of 

the  spin  dependent  terms)  tend  to  diminish  to  zero  when  the  system  Is  In  a 
disordered  state  and  tends  to  grow  when  the  system  Is  In  an  ordered  state. 

For  a  given  Initial  set  of  €%*  therefore  adjust  ji  to  obtain 

the  order-disorder  transition  temperature  from  the  behavior  of  p  and  the 
coupling  constants  after  many  RG  transformations. 

After  the  th  RG  transformation,  the  free  energy  |5PG  of  Eq(17)  may 
be  approximated  by 

,(«<) 


p  (  wax  ^)  ) 


(21) 


when  .<  Is  sufficiently  large.  Eq(21)  contains  the  contribution  to  the  free 
energy  from  the  spin-independent  constant  terms  in  each  as  well  as  an 

( 

approximation  for  the  contribution  from  the  spin-dependent  part  of  We  can 

also  calculate  the  first  and  the  second  derivative  of  ^PG  with  respect  to 
temperature  and  thus  obtain  the  internal  energy  and  specific  heat  of  the  system. 

Such  calculations  involve  evaluating  the  product  of  a  series  of  matrices 

'  dV'  I  differentiating  numerically  and  will  be 

detailed  in  another  paper. 

In  discussing  the  MKVM  we  have  neglected  the  coupling  between  the  spins 
in  the  R  and  R*  sublattices  completely  up  to  now,  i.e.  we  have  let  0  in  the 

last  term  in  Eq(7).  As  argued  in  Section  2,  this  term  is  unimportant  near  T^. 
Therefore  we  may  take  its  effect  into  account  approximately.  Fq(ll)  may  be  rewritten 


^  (O' ]  (22) 

■*  ft' 


The  coupling  of  spins  the  fir'-  bracket  is  shown  in  fig.  3. 
interaction  energy  V  (Tl'  n  eai.h  cell,  which  has  four  spins 

K  *N 


In  addition  to  the 
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In  the  corners,  there  is  one  additional  spin,  say  ,  in  the  center, 
which  interacts  with  each  corner  spins  with  coupling  constant  6,  /g'  . 

These  central  spins  actually  interact  with  each  other  via  the  first  term 
in  the  second  bracket  of  Eq(22).  In  our  approximation,  we  neglect  this 
interaction  between  the  spina  in  the  first  and  the  second  brackets. 

We  consider  the  effect  of  the  central  spin  in  the  first  bracket  as  being 
only  to  produce  a  certain  effective  coupling  among  the  corner  spins.  In 
this  approximation,  we  still  have  two  decoupled  sublattices.  We  carry  out 
the  configuration  sum  for  the  central  spin  first  and  obtain  the  effective 
cell  potential  for  the  corner  spins: 

<  I  h  f  €  ,<p(  *p(  -P  f ^  ^ 

We  now  obtain  effective  initial  values  of  ^  the  four-spin  coupling 

and  the  spin-independent  constant  term  by  comparing  Eq(23)  and  Eq(8).  The  RG 
calculation  procedure  can  then  be  carried  through  as  in  the  =  0  case 

described  above. 

4.  Results  for  0/W(llQ) 

We  employed  the  method  introduced  in  Section  3  to  calculate  the  order- 

disorder  transition  temperature  for  several  sets  of  ^  values 

used  in  the  Monte  Carlo  calculations  of  E.D.  Williams  et  al.  The  results 

are  shown  in  table  1.  Columns  1,  2,  and  3  list  the  interaction  energies  arranged 

according  to  the  ratio  -l«t  /fi,  which  is  given  in  Column  4.  Column  5  gives  the 

-> 

order-disorder  transition  temperature  T^  obtained  by  Williams  ct.al.  using  the 
Monte-Carlo  method.  Column  6  given  the  value  of  T  for  =0  in  the  free  fermion 
approximation  of  Fan  and  Wu.  Our  results  arc  nho\^n  in  Column  7  where  the  effect 
of  Is  neglected  and  column  8  where  the  effect  of  is  included  as  described 
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/T^  (rue) 

in  Section  3.  The  ratios  shovm  in 

Columns  9  and  10,  respectively.  From  table  1  we  find  that  for 

~  ^•’7  the  difference  between  T  (ter  0)  and  T  (MC)  is  always  less 

c  T  c 

than  8%  and  for  —  ^  oM'  the  difference  between  T(4-, 0)  and  T  (MC) 

'  c 

is  always  less  than  10%.  Thus  our  method  quite  accurately  predicts  the  transition 
temperature  at  half  coverage  for  a  certain  range  of  •  This  range  is  close 

to  physical  reality  because  1b  a  longer  distance  interaction  than 

and  thus  will  be  weaker  in  general.  It  is  also  clear  from  Table  I  that 
while  our  approximate  handling  of  increases  the  error  in  T^  somewhat,  the  major 
part  of  it  is  due  to  the  underlying  treatment  of  the  next  nearest  neighbor  Ising 
model. 

We  have  also  calculated  the  average  internal  energy  and  specific  heat 

of  the  system  for  the  set  of  data  ^,=  -2.0  KCal/mole,  ^2' =  1.74  K  cal/mole 

and  =  -0.9  K  cal/mole  without  taking  the  effect  of  last  term  in  Eq(7)  into 

account.  The  results  are  shown  in  Fig.  4.  Comparing  Fig.  4  with  Fig.  7  of 

ref.C'Zj,  which  was  obtained  by  Monte-Carlo  method,  it  is  easy  to  see  that 

our  internal  energy  is  very  close  to  that  obtained  by  the  Monte-Carlo  method. 

However  the  percentage  difference  between  our  specific  heat  and  that  obtained 
li«? 

by  Monte-Carlo  method  is  relatively  large  for  temperatures  far  from  T^. 

Our  specific  heat  is  larger  than  that  obtained  by  the  Monte  Carlo  method  for 


5 •  Critical  Properties 

The  critical  properties  of  the  0/W(110)  system  are  of  considerable 

interest.  This  is  hecaufic  .an  argument  based  on  the  s>Timetry  of  the  (2x1) 

ordered  state  puts  the  order-disorder  transition  in  the  universality  class  of  the 
\7  (  iaJ.I-L  (.‘W.l  O.'J  iti- 

,^eight  vertex  nodel'."  This, model,  for  which  an  exact  solution  exists,  exhibits 

h  K 

continuously  varying  critical  order.  That  is,  the  transition  is  characterised  by 
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crltlcal  exponents  that  vary  continuously  as  a  function  of  certain  Interaction 
parameters.  In  a  renormalization  group  analysis,  the  mechanism  for  this  Involves 
a  line  of  fixed  points.  ' 

Now  for  0/W(110)  there  has  been  no  measurement  of  critical  properties 
to  date.  Further,  the  argument* 'that  they  should  be  In  the  eight  vertex 
universality  class  Is  based  on  the  symmetry  of  the  ordered  state  only,  and 
involved  some  assumptions  that  have  not  been  completely  justified.  So  the 
question  of  the  critical  properties  of  this  system,  or  Eq(l) ,  is  of  considerable 
interest. 

We  present  several  facts  that  indicate  that  the  order-disorder  transition 

of  Eq(l)  is  Ising-llke  at  one  half  coverage.  This  would  mean  that  this  coverage 

il 

corresponds  to  the  decoupling  point  (four  spin  coupling  K  =  0)  of  the  eight 
vertex  model.  Assuming  the  group  theoretic  arguments  mapping  0/W(110)  onto  the 
eight  vertex  model  are  correct,  this  would  imply  that  the  exponents  along  the 
rest  of  the  order-disorder  phase  boundary  vary  continuously  away  from  the  Ising 

A  Z 

values.  (Note  especially  that  at  1/4  is  not  independent  of  ). 

First  we  note,  in  Eq(l),  that  if  ^|=  0  the  model  reduces  to  a  simple 

antiferromagnetic  Ising  model  with  next-nearest  neighbor  interactions  (on  each 

Sublattice).  Such  a  model  has  not  been  solved  exactly,  but  there  Is  considerable 

evidence  that  it  is  in  the  universality  class  of  the  simple  Ising  model.  Now  as 

Z 

pointed  out  previously,  Monte  Carlo  results  show  that  the  transition  temperature 

I, 

at  one-half  coverage,  is  independent  of  6^  over  a  considerable  range  of 

values.  We  have  shown  above  that  this  remains  true  oven  if  one  considers 

(see  Table  I),  at  least  within  expected  error.  iiThite  it  is  possible  for  t  j  to 

1, 

affect  the  critical  exponents  without  changing  T  it  is  quite  unlikely,  since 
this  would  require  a  line  of  fixed  points  parallel  to  the  over  a  very 

large  region.  Die  oighr  vertex  fixed  line  is  such  that  T^  depends  on  the  "distance” 
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from  the  decoupling  point.  This  picture  Is  also  consistent  with  the 
demonstration.  In  Section  2,  that  the  6,  term  Is  Irrelevant  at  the 
6,  •=  0  fixed  point.  ' 

A  second  bit  of  evidence  Is  that  we  are  able  to  get  good  results  for 
with  a  renormalization  group  treatment  that  has  only  one  (Islng-llke) 
fixed  point.  Note  that  we  do  Include  6-^  in  our  calculation,  but  its  only 
effect  is  to  renormalize  the  other  coupling  constants. 

Z 

Finally,  we  recall  the  argument  of  Williams  et.al.  mentioned  in  section  2. 
Consider  a  perfectly  ordered  (2x1)  state  (at  6  =  h) •  Moving  an  adatom  to  the 
most  favorable  disordered  site  raises  Its  energy  by  At  ~  or 

AC -Itj  f  >  depending  on  the  orientation  of  the  ordered  state.  This  At- 

h 

is  independent  of  ,  which  Is  again  consistent  with  the  Independence  of 

and  an  Ising  type  transition  since  this  model  is  Ising  like  for  =0. 

II 

Now  consider  a  spin  representation  of  the  eight  vertex  model.  In  the 
perfectly  ordered  state,  the  tiro-spin  excitation  that  preserves  the  number  of 
up  spins  and  has  smallest  At"  in  general,  involves  both  sublattices.  However 
when  the  four-spin  coupling  K  =  0,  this  is  no  longer  true.  Then  one  can 

^ ( (  II  t 

^Gont-inue  the  excitation  to  a  single  sublattice  of  the  model,  which  has  simple 
Ising  coupling. 

6.  Future  Rc.search 

It  is  clear  from  the  results  reported  here  that  the  renormalization 
group  technique  is  a  u.seful  method  for  extracting  AA  interaction  values  from 
experimental  .surface  phase  diagrams  with  reasonable  accuracy.  To  complete 
this  program,  however,  it  is  necessary  to  generalize  the  method  to  coverages 
Cr  *5.  We  have  done  some  preliminary  work  and  find  reasonably  good  agreement 
for  at  0  =  '4  for  the  AA  interaction  parameters  ^  i' 

kcal/mole.  However  it  is  clear  that  an  accurate  treatment  for  6*  j*  -j  requires 
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generalizing  the  RG  tranoformatlon  in  Eq(19)  to  Include  a  term 

in  the  exponent,  l.e.  an  effective  magnetic  field  (or  staggered  field) 

to 

acting  on  the  transformed  spin  variables.  Finding  values  of  and  p 
that  give  accurate  results  will  require  some  further  study. 
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Figure  Captions 


Fig.  1. 

Fig.  2. 
Fig.  3. 
Fig.  4. 


Table  I 


Absorption  sites  and  adatom-adatom  in;teractlons  for 
0/W(110) .  Both  •  and  x  denote  adsorption  sites.  The  former 
comprise  the  R  sublattice  and  the  latter  the  R'  sublattice. 

Location  of  four  spins  on  a  typical  R  cell  of  the  lattice. 

Spin  couplings  for  the  first  bracket  of  £q(21). 

Of 

Internal  energy  E  and  specific  heat  for  ~  2.0,  “  1-74 

and  =  -  0.9  kcal/mole  calculated  by  the  MKVM  method. 


Table  Caption 


.  Columns  1,2,3:  adatom-adatom  interaction  energies 

Column  4:  -  c, /e  . 

4  I 

Column  5:  transition  temperature  T^  as  determined  by  Monte 
Carlo  calculation^ 

‘J 

Column  6:  T^,  from  the  free-fermion  approximation  (for  the  case). 

Columns  8  and  9:  Tq  from  the  present  work  with  £^^=0  or  respectively 

Columns  9  and  10:  ratios  of  the  present  values  of  T^,  to  the 
Carlo  results^ 


Monte 
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